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I. INTRODUCTION 



During the last years there has been a tremendous research area on low dimensional 
semiconductor systems such as quantum wells, quantum wires and quantum dots (QDs).->2 
Low dimensional semiconductor systems display various physical properties, different from 
higher dimensional systems for instance, the discrete energy structure, the absorption spec- 
trum (both intra- and inter-band) are expected here to be a series of discrete lines.- 

Recent improvements in the nanostructure technology have made possible to prepare zero 
dimensional semiconductor nanostructures, called quantum dots, by using different meth- 
ods such as molecular-beam epitaxy (MBE). 4 As it is well known, when dimensions of the 
structure are diminished, the quantum mechanical effects become more important and hence 
these effects improve the performance of many electronic and optoelectronic devices. The 
attracting features of QDs are that their many properties (energy levels, wave functions, den- 
sity of states, etc.) can be controlled by the shape, size and composition of the structure. 
These features are similar to the properties of real atoms^ and therefore QDs are named as 
artificial atoms. A great number of studies on quantum dot systems such as GaAs/AlGaAs, 
InGaAs/GaAsS, InAs/InP-^, and InAs/In(Ga,As)2 have been reported both experimentally 
and theoretically for a wide range of dot sizes and shapes. The studies on QDs open a new 
field in both basic physics and chemistry, and present a wide range of potential applica- 
tions for optoelectronic devices such as optical and electro-optic modulators 1 ^ 1 -, inter-band 
lasers 12 , inter-subband long wavelength detectors.— 

A hydrogenic impurity problem in a QD is a very useful model in understanding of 
many electronic and optical properties of these kind of structures. A donor impurity atom 
has a one more electron than that required to make a chemical bond with neighboring 
atoms in a semiconductor materials. If a donor impurity has an additional one electron, 
such impurity is known as hydrogenic one because of that it is very similar to a hydrogen 
atom. 14 The presence of an impurity in a QD changes its effective potential and this affects 
both energy spectrum and optical transitions because of the Coulomb interaction between 
the electron and the impurity and is very important in low dimensional semiconductor 
physics.— Also, some device properties, such as transition energy, can be tuned by the 
impurity. Therefore, the effects of impurities on the optical and electronic properties of 
semiconductor heterostructures draw attention most of the authors and so a number of 



2 



theoretical and experimental studies have been reported.— — 

As is well known, when an electron is stimulated by a photon, the electron makes a 
vertical transition from an initial level to a final one. If the levels are in the same band 
such as conduction band, this kind of transition is called as inter-subband or intra-band 
transition. On the other hand, if the transitions occur in between valance and conduction 
bands, these transitions are described as inter-band ones. Similar kinds of transitions take 
place in QDs too. Today, the most of optoelectronic devices fabricated from QDs have 
been built on such transitions. Hence, investigation of the optical absorption properties of 
QDs are very important in condensed matter and applied physics.— - — Li and Xia^ 5 -, and Li 
et aL 26 have calculated both inter-band and intra-band optical transitions for InAs/GaAs 
coupled quantum dots. Han et al.— have calculated the inter-subband optical transition of 
QDs in a quantum well in the framework of the effective mass envelope function theory. 
Yilmaz and Safak 28 , have examined the oscillator strengths of inter-subband transitions for 
the electron in CdS/Si02 QD with an on-center donor impurity. Bassani and Buczkc 29 
have showed that how to determine the resonant transition and their lifetimes and how 
to calculate the oscillator strengths for inter-sublevel optical transitions from the ground 
state. De Souso et al.— have studied the influence of the Stark effect on the intra-band 
transitions due to its potential application in infrared detection and emission. The band 
structure and inter-sublevel hole transitions in valance band of SiGe/Si QDs have been 
investigated theoretically by Lin and Singh 3 -^, using eight-band kp model. Gondar and 
Comas^ 2 - have obtained the selection rules for the optical transitions and they reported the 
oscillator strengths for dip ole- allowed transitions in a QD with semi-spherical geometry. 
Jiang et al.— have performed theoretical studies on inter-sublevel transitions in two QD 
systems. 

The recent developments in fabrication technology make possible the production of multi- 
shell quantum heterostructures.— The experimental studies carried out on these structures 
are mostly fulfilled based on inter-band exciton transitions. 35 However, although there are 
some theoretical studies reported on MSQDs in the literature, the studies are related to the 
electronic properties of these structures such as energy levels and impurity binding energies.- 3 ^ 
Recently, Sahin et al~^ has showed that the electronic shell structure of MSQD is changing 
and reordering with shell thicknesses. They have demonstrated that the reordering in the 
electron-conduction band shell structure can be completely different from the hole-valance 
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FIG. 1. Schematic representation of a core/shell/well/shell QD and its potential profile. 

band shell structures. In our previous work 37 , we have performed a detailed investigation 
of the electronic properties of a core/shell/ well/shell multi-layered spherical quantum dot 
for both ground and excited states. To the best of our knowledge, the optical properties of 
MSQD have not been investigated yet. The controllability of the layer(s) thicknesses may 
provide different advantages in the production of new generation opto-electronic devices 
working based on inter-sublevel optical transitions. 

The main goal of this study is to research the optical properties such as inter-sublevel ab- 
sorption and the oscillator strength between ground and excited states of a MSQD structure. 
These properties have been studied for cases with and without an on-center hydrogenic donor 
impurity. The calculations have been done for various core radii, barrier thicknesses and 
well widths for same potential confinement. All results have been presented and discussed 
comparatively for both cases with (Z = 1) and without the impurity (Z = 0). 

The rest of the paper is organized as follow: In the next section, we introduce the consid- 
ered model and its theory. In section IHIl the results of the calculations and their probable 
physical reasons are discussed. In the last section, a brief conclusion is presented. 

II. MODEL AND CALCULATIONS 

In our study, we deal with a system which includes an electron confined in a MSQD. 
The considering system consists of two GaAs quantum dot one within the other and each 
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surrounded by AlGaAs layers. The first AlGaAs layer generates a finite potential and it 
allows the electron to tunnel between the core and well regions. The other AlGaAs layer 
isolates the whole system from outside. The schematic representation of MSQD and its 
potential profile are plotted in Fig. [TJ In the considering structure, the core radius is Ri, 
the shell thickness is T s = R 2 — R\ and the well width is T w = R 3 — R 2 . In the effective mass 
approximation for a spherically symmetric quantum dot, the single particle Schrodinger 
equation is written as 



— V ( 1 V 1 - + + H2 + v(r) 

2 r \m* (r) r J n{r)r 2m*{r)r 2 



Rnf (r) = E n tR n t (r) 



Here, h is reduced Planck constant, m*(r) is the position-dependent electron's effective 
mass, Z is the charge of the impurity, k(t) is the position dependent dielectric constant, 
i is the angular momentum quantum number, V(r) is the finite confining potential, E nt g 
is the electron energy eigenvalue for specified principle quantum number (n) and angular 
momentum quantum number (£), R n /(r) is the radial wave function of the electron. It should 
be noted that the Z = and Z = 1 correspond to a case without and with a hydrogenic 
donor impurity, respectively. The mathematical expression of the confining potential is 



' 0, r < R 1 and R 2 < r < R 3 



V(r) 



(2) 

Vb, R\ <r < R 2 and r > R 3 

In order to determine the single particle energy levels and corresponding wave functions, 
Eq.(JT]) is solved full numerically by the shooting method. As is well known, this technique 
converts an eigenvalue problem to an initial-value problem. For this purpose, Hamiltonian 
operator is discretized on a uniform radial mesh in ID using the finite differences, then 
Eq.([Tl) can be reduced to an initial- value equation by means of 
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where i is the index of mesh points, h is the distance between two mesh points and it is 



chosen as 0.005. The details of this method can be found in Ref. 



14 



5 



The photon absorption process can be described as an optical transition that takes place 
from an initial state to a final one with assisted by a photon. The optical absorption 
calculations for the inter-sublevel transitions are based on Fermi's golden rule derived from 
time dependent perturbation theory. The inter-sublevel optical absorption coefficient is 
given as ■ 

a (hw) = 167r2(3 ^ sNif hu\z fi \ 2 5(E f -Et- tkj), (4) 

"T 'con 



where n r is the refractive index of the semiconductor and this value is 3,15 for GaAs , V con is 
volume of the confinement regions, fruj is the incident photon energy, ftps is the fine structure 
constant and its value is 1/137, Nif — Ni — Nf is the difference of the number of electron 
between the initial and final state, and Ei and Ef are energy eigenvalues of the initial and 
final state, respectively. Zif is the dipole matrix element between the initial and final states 
and the S is broadening parameter. 

In spherical quantum heterostructures, the selection rule, A£ = ±1, determines the final 
state of the electron after the absorption. Hence, the Is state is taken as the initial state 
and the lp level is considered as the final state. The total wave function of these levels are 
determined from multiplication of the radial wave functions with the spherical harmonics 
(i.e. Rrd(r)Y tm (0,<t>)). 

The dipole matrix element for transitions between Is — lp levels is 



\ Z f: 



2 

,2 1 ' ^ 
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R f (r)r 3 Ri(r)dr 



o 



(5) 



Ri(r) and Rf(r) are radial wave functions of the initial and final states, respectively. The 
factor | comes from the integration of the spherical harmonics. In addition, the 5 function 
in absorption equation is replaced by a narrow Lorentzian by means of 

5(E f -E i -hu) = — — 5 5-, (6) 

1 f ' vr [(Hlj - (E f - E t )) 2 + (hT) 2 } K ' 

where the term of hT is maximum width in semi height of the Lorentzian and its value is 

taken 6.4 meV in this study. 

Other quantity in the study of optical properties is the determining of oscillator strengths. 

The oscillator strength gives the information about magnitude of the absorption. That is, 



the amount of the oscillator strength is directly proportional to the absorption coefficient. 
For Is — lp transitions the oscillator strength is 



2m* 



R f (r)r 3 Ri(r)dr 



(7) 



III. RESULTS AND DISCUSSION 

The atomic units have been used throughout the calculations, where Planck constant 
h = 1, the electronic charge e = 1 and the bare electron mass m = 1. Effective Bohr 
radius is do — fOOA and effective Rydberg energy is R y ~ 5.25 meV. The material param- 
eters have been taken as m Gaj 4 s =0.067m , m^ Ga ^ s =0.088m , H=228 meV, K Gaj 4 s =13.18, 
KAiGaAs= 12. 8. Also the effective masses of electrons inside GaAs and AlGaAs are ml and 
an d the dielectric constants are n\ and k 2 , respectively. The position-dependent effective 
mass and the dielectric constant may be defined as follows^ 



m*(r) 



r < R 1 and R 2 < r < R 3 



^i, Rx<r <R 2 and r > R 3 



1, r < i?i and R 2 < r < R 3 



(8) 



^, R l < r < R 2 and r > R 3 

In Eq. HI all terms except Nif, Jtlu, V con and Lorentzian are constants and these constants 
have no effect on the determining of absorption shape. In the resonant absorption case which 
occurs when the photon energy becomes equal to the energy difference between the levels, 
the Lorentzian is also a constant, 1/itKT. As the optical absorption coefficient is directly 
proportional with the dipole matrix element, it is inversely proportional with the volume 
of the confinement regions. Hence, the overlapping, volume of the confinement regions 
and the energy levels become more effective on the absorption coefficient. Therefore, the 
probable physical reasons in the absorption properties will be discussed as dependent on 
these parameters. Similarly, the oscillator strength also depends on the energy difference 
between considered levels and the overlapping of the wave functions of these levels. 




FIG. 2. (Color online) Variation of the absorption spectra as a function of incident photon energy 
for different core radii in cases of Z = and Z = 1. The shell thickness is T s = 0.2 ao, the well 
width is T w = 0.4 ao- 



It should be noted that the overlapping of the radial wave functions is calculated by 
means of 



A = / R f ( r )Ri(r)r 2 dr. 



(9) 
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FIG. 3. (Color online) Ground and excited states probability distributions as a function of the r 
for different core radii. The top panel for Z = and the bottom panel for Z = 1 case. The energies 
of levels are showed with dashed lines. 

A. Effect of the core radius on the optical properties 

Figure [2] shows the absorption coefficient variation with the incident photon energy for 
Z = and Z = 1 cases for different core radii and constant shell thickness (T s ) and well 
width (T w ). As it is known, when an electron is stimulated by a photon, the electron absorbs 
the photon and moves up in the vertical direction to higher energy levels. If the incident light 
energy equals to the energy difference between the levels, this energy is called as resonance 
transition energy which corresponds to the peak energy of the absorption coefficient. In 
Fig. [21 the peak energies of the absorption coefficient for Z = 1 case have higher values 
(blue shift) than those of Z = case at all core radii. This situation results from the 
attractive Coulomb potential of the donor impurity. Because, in case of Z = 1, the impurity 
pulls ground (Is) level down more than excited (lp) level and hence the energy difference 
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between Is and lp levels becomes larger comparing with Z = case as seen evidently from 
energy levels of Is and lp states in Figj3l On the other hand, although the absorption 
coefficient values are almost same for Z = and Z — 1 cases at R\ = 0.4 and 0.8 ao, it is 
larger for Z = than in case of Z — 1 for Ri = 0.6 a . This can be explained by means 
of the overlapping of the wave functions given in Fig. [3] for Z = and Z = 1 cases. That 
is, while the overlapping of the wave functions are almost same in both Z = and Z = 1 
cases for Ri = 0.4 and 0.8 ao as seen from right and left panels of Fig. [3j the overlapping is 
stronger in case of Z = than that for Z = 1 case because the finding probability of ground 
and excited levels are confined in different region in latter case by comparison with Z = 
case as seen from middle panel of Fig. [3J 

The variation of peak values of the absorption coefficient with the core radius may provide 
some important information in device applications. For this purpose, the absorption peak 
values are plotted as a function of R\ for Z = and Z = 1 cases in Fig. HJ The structure 
parameters are given on the figure. According to the figure, the absorption coefficient 
presents the downward tendency in general for both cases except that sudden decreasing 
and increasing. This sudden decreasing and increasing in the absorption coefficient with Ri 
can be explained as depending on overlapping of the wave functions. When the R\ = 0.3 ao 
ground and excited states wave functions are localized in the well region. Therefore, the 
overlapping is strong and hence the dipole matrix element has a big value. In this situation, 
the optical transition takes place in the well region. With the increase of the core radius, 
the overlapping decreases and so the absorption coefficient decreases too. In this case, while 
the wave function of ground state is localized in the core region, the excited states wave 
function is confined in the well one. In large core radii (_Ri > 0.9 a ), the wave functions 
of ground and excited states are localized in the core region completely and therefore the 
absorption coefficient becomes large for the cases of Z = and Z = 1. The absorption takes 
place in the core region anymore. On the other hand, although the overlapping increases, 
the decrease in the absorption coefficient can be explained by the increase of the confinement 
region volume which results from the increase of R\. 

The oscillator strength variation of the transition as a function of the core radius is given 
in Fig. |5] for Z = and Z = 1 cases. The inset figure shows the overlaps of the wave 
functions calculated from Eq. |9j From the figure, we see that the oscillator strength does 
not show a considerable change in small core radii. After that, it exhibits a sudden fall 
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FIG. 4. (Color online) Absorption coefficient peak values as depending on core radius for Z = 
and Z = 1 cases. The shell thickness is T s = 0.4 ao, and well width is T w = 0.4 ao- 

and rise tendency similar as the absorption coefficient. It is remains fixed about 1.0 after 
Ri > l fl o with further increasing of R\. The oscillator strength is strongly dependent on 
both energy difference between the levels and overlapping of the wave functions as it is 
mentioned before. Although the energy difference is little in small values of R\, overlapping 
of the wave functions is strong as seen from the inset and thus the oscillator strength is 
also higher. The overlapping exhibits a decrease tendency rapidly with increasing Ri. This 
decreasing situation is faster in Z = 1 case. In contrast, when R\ > 0.6 ao, the increase in 
overlap of Z = 1 case becomes more quickly than that of Z = 0. In addition to this, the 
energy difference gets some more in Z = 1 case as seen from Fig. HI In parallel with these, 
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FIG. 5. (Color online) The variation of oscillator strength as a function of core radius for Z = 
and Z = 1 cases. The structure parameters are same with the previous figure. The inset figure 
shows the overlap of the wave functions. 

the oscillator strengths of the impurity are larger after certain Ri value. 

B. Effect of the shell thickness on the optical properties 

Fig. [6] shows the variation of absorption coefficient spectra as depending on the photon 
energy for different shell thicknesses and the other parameters of considering structure are 
Ri = 0.6 ao and T w = 0.4 a^. The figure is comparatively plotted for Z = and Z = 1 
cases. As seen from the figure, the absorption coefficients of Z = case are higher than that 
of Z — 1 case for all shell thicknesses. When the shell thickness increases, the absorption 
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FIG. 6. (Color online) Absorption coefficient as a function of photon energy at different shell 
thicknesses for Z = and Z = 1 cases. Other structure parameters are specified on the figure. 

coefficient decreases very rapidly for both Z = and Z = 1 cases. Therefore, the absorption 
coefficients are multiplied by 2 and 10 for T s = 0.4 and T s = 0.6 ciq, respectively. In addition, 
although there are no significant changes in the absorption peak energies (i.e. transition 
energy) for Z = 1 case, these changes are notable and downward tendency with increasing 
shell thickness in case of Z = 0. The probability distributions of dealing structure for Z = 
and Z = 1, which are top and bottom panel of Fig[7J respectively, are given for T s = 0.2, 
0.4, and 0.6 a>o for certain core radius and well width. As seen from the figure, ground and 
excited states wave functions become localized in different region with the increase of the T s . 
For T s = 0.2 ao, while a part of ground state finding probability is in the well region in case 
of Z = 0, because of attractive Coulomb potential of the impurity it is almost completely 
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FIG. 7. (Color online) Ground and excited states probability distributions as a function of r for 
different shell thicknesses in case of Z = (top panel) and for Z = 1 (bottom panel). 



in core region for Z = 1 case. On the other hand, the finding probability of excited states 
are localized in the well region for both Z = and Z — 1 cases. This is because the 
probability of tunneling, between core and well regions, decreases with the increase of the 
shell thickness. This process makes weak the overlapping of the wave functions of ground 
and excited states. Also, when we look at the top panel of Fig. 0, we observe that the 
energy difference between the Is and lp states decreases with the increase of shell thickness 
and hence, the peak energies of the absorption coefficient decreases (red shift) as seen in Fig. 
EJ When we compare the absence of impurity with the existence of one, in Z = 1 case, the 
impurity pulls the Is energy states down, thus the energy difference between ground and 
excited levels is larger than that of Z = case. On the other hand, the energy levels remain 
almost fix with increasing of the shell thickness for Z = 1 case and hence the absorption 
peak energies become steady. 
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FIG. 8. (Color online) The variation of peak values of the absorption coefficients as depending on 
T s for Z = and Z = 1. The core radius is R\ = 0.6 ao, and the well width is T w = 0.4 ao- 

Fig. [8] shows the variation of the maximum values of the absorption coefficient with the 
shell thickness, for Z = and Z — 1 cases. The parameters of the structure are given on 
the figure. As is expected, the peak values of absorption coefficient exhibit a downward 
tendency with the increase of T s . Although there is no change in the confinement region 
volume, the decrease of the overlapping lead to this downward tendency. Whereas, in small 
shell thicknesses, the difference of absorption coefficient for with and without the impurity 
is very small, this difference rises gradually and remains almost constant in large shell 
thicknesses. Moreover, the values of absorption coefficient for Z = 1 are lower than that of 
Z = 0. In this case, the decrease of the overlapping becomes more effective than the increase 
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FIG. 9. (Color online) Oscillator strength variation as depending on T s for Z = and Z = 1. The 
core radius is R\ = 0.6 do, and the well width is T w = 0.4 ao- The inset shows the overlap of the 
wave functions. 

of the energy difference between the states for Z = 1. 

The oscillator strength of the structure as a function of T s is given in Fig. [9] for the 
same parameters with the previous figure. The inset figure shows the variation of the 
overlap integral with the shell thickness. As seen from Fig. [HJ the oscillator strength has a 
maximum value at the beginning and reduces rapidly with the increasing of T s , and it goes 
towards zero with further increasing of the T w for both Z = and Z = 1 cases. When we 
look at the figure, it is observed that the oscillator strength for Z = 1 smaller than that for 
Z = until a certain value. This case stems from the impurity atom. Because, the impurity 
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FIG. 10. (Color online) The variation of absorption coefficient as a function of the photon energy 
for different well widths in cases of Z = and Z = 1. The core radius is R\ = 0.6 ao, and the shell 
thickness is T s = 0.2 ao- 



pulls the ground state down and so, it localized in the core region while the excited state 
localized in the well one. Thus, as seen from the inset, the overlapping of the wave functions 
and as depending on this, the oscillator strength decreases. Also, the increase in T s causes 
the decrease in tunneling between the core and well regions. Therefore, these regions do not 
feel each other. 
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FIG. 11. (Color online) Ground and excited states probability distributions as a function of r for 
various well widths in cases of Z = (top panel) and Z = 1 (bottom panel). The other structure 
parameters are specified on the figures. 



C. Effect of the well width on the optical properties 



The variations of optical absorption coefficient with the well widths may be important for 
different device applications of MSQDs. For this reason, we have calculated the absorption 
coefficients depending on the well widths. Figure [10] shows the absorption coefficients as 
a function of incident photon energy in case of Z = and Z — 1 for three different well 
widths, i.e. T w = 0.2, 0.4 and 0.6 ao- The core radius and the shell thickness are chosen 
as Ri = 0.6 a and T s = 0.2 a , respectively. When we look at the figure, we observe that 
although the absorption coefficients of Z = bit larger than that of Z = 1, 

their peak energies of Z = 1 case are greater (blue shift) than that of Z = for given well 
widths. Also, it can be seen that the absorption coefficients and their peak energies are 
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higher for lower T w and reduce with increase of the T w for both Z = and Z = 1. This can 
be explained as follows with the assistance of Fig. [TT] since the well width is narrow, the 
probability densities of the levels are higher in the core region and hence, the overlapping of 
the wave functions increases relatively with respect to larger well widths (i.e. T w > 0.4 ao) 
and it leads to the absorption occurs in the core region as seen from Fig. [TT] (a) and (d). In 
this case, the energy difference between the states is also bigger. In addition, the volume of 
the confinement region is small due to the narrow well width. As a results of these, when 
T w = 0.2 do, both the absorption coefficients and their peak energies become larger. When 
T w = 0.4 a , the absorption coefficients and their peak energies reduce. This is because the 
energy level of excited state decreases considerably owing to the tunneling of excited state 
wave functions to the well region for both Z = and Z = 1 cases as seen in Fig. [TT] (b) and 
(e). In case of Z = 1, because of attractive Coulomb potential of the impurity, ground state 
wave functions almost completely confines in the core region as excited state wave function 
confines in the well one. This results in small overlapping of the wave functions and so small 
absorption coefficient as compare with case of Z — 0. When T w = 0.6 ao, the probability 
density of excited states are confined completely in the well region in both Z = and Z = 1 
cases as seen from Fig. [TTJ (c) and (f). On the other hand, ground state density of Z = is 
localized almost completely in the well region as that of Z = 1 case is shared between the 
core and well regions. When T w > 0.6 a , the absorption takes place in the well region. 

Fig. [12] shows the variation of absorption coefficient peak values as a function of the T w 
in cases of Z = and Z = 1 for R\ = 0.6 ao and T s = 0.4 ao- As expected, the absorption 
coefficient exhibits a reduce tendency with the increase of the T w in general. When T w is 
about 0.4 ao for Z = 0, the absorption coefficient becomes minimum and it increases again 
with further increasing of T w . In Z = 1 case, the absorption coefficient reaches its minimum 
value when T w = 0.5 a . These decreases can be explained by means of that the overlapping 
of the wave functions is weak as seen from the inset given in Fig. [13] In Fig. [12] it is 
observed that the absorption peak values in between T w = 0.3 and 0.7 ao are smaller for 
Z = 1 case by comparison with Z = case. Because, the attractive Coulomb potential 
of the impurity reduces the overlapping of the wave functions. When T w = 0.7 ao, the 
absorption coefficients of both cases have same values and exhibit same treatment. That is, 
both the core region and the impurity lose of their own effects on the considered quantum 
states anymore and the optical transitions take place in the well region. 
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FIG. 12. (Color online) The variation of absorption coefficient peak values as depending on T w for 
Z = and Z = 1. The core radius is R\ = 0.6 ao and the shell thickness is T s = 0.4 ao- 

In Fig. Q3J for the structure with the i?i = 0.6 ao and T s = 0.2 ao, the variation of the 
oscillator strengths are given for cases with and without the impurity depending on T w . It 
is seen that, in the beginning, the oscillator strengths have maximum values and then, they 
exhibit a decrease tendency with increasing of the T w and reach minimum values. In the 
latter case, ground state wave function is localized in the core region, while excited state 
wave function is confined in the well one. Afterwards, they increase with further increasing 
of the T w and decrease again smoothly. As similar to the previous figure, when almost 
T w < 0.7 a , the oscillator strength in the presence of the impurity has smaller value due to 
the negative effect of the impurity coulomb potential on the overlapping. After this certain 
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FIG. 13. (Color online) The oscillator strength variation as depending on T w for Z = and Z = 1 
cases. The core radius i?i = 0.6 ao, and the shell thickness is T s = 0.2 ao- The inset shows the 
overlap of the wave functions. 



value of T WJ the oscillator strength is to become higher for Z = 1 than that of Z = case, 
although the overlapping is equally stronger in both observed from the inset of Fie;. 

fl~31 The reason of this alteration is because the difference between the energy levels is larger 
in case of Z = 1 owing to attractive coulomb potential of the impurity on ground state, 
especially. We conclude from here, the energy difference becomes more predominant on the 
oscillator strength than the overlapping. 
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IV. CONCLUSION 



In this study, we have calculated the absorption coefficient and oscillator strength for 
the inter-sublevel transitions in a MSQD heterostructure. Here, Is and lp levels are taken 
into consideration as ground and first excited states, respectively. The results have been 
evaluated depending on core radii (Ri), shell thicknesses (T s ), and well widths (T w ) and 
we have discussed the probable physical reasons behind the results. We have observed that 
the optical properties are drastically dependent on the layer thicknesses of the structure. 
We have concluded that many of the optical properties can be controlled by tuning the 
layer thicknesses for different device applications. To our knowledge, this study is first 
investigation of the detailed optical properties of a MSQD. We believe that this study will 
be rather useful and will contribute to the understanding of the optical properties of MSQDs. 
We hope that this study will also stimulate both experimental and theoretical investigations 
of the optical properties of MSQDs. 
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